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Abstract. The 16ms X-ray pulsar PSR J0537-6910 in 
SNR N157B is the fastest known isolated (non-recycled) 
pulsar and, with a rotational energy loss E ~ 4.8 x 10 38 
erg s _1 , is the most energetic (together with the Crab). 
Here we report the results of optical observations of the 
field, recently performed with the SUSI2 camera of the 
NTT. Few objects are observed inside the ~ 3" X-ray 
error circle but none of them can be convincingly as- 
sociated to the pulsar, which appears undetected down 
to V ~ 23.4. With a corresponding optical luminosity 
L opt < 1.3 x 10 33 erg s _1 , PSR J0537— 6910 is, at best, 
comparable to the other very young pulsars Crab and 
PSRB0540-69. 
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1. Introduction 

PSR J0537— 6910 is a young, fast, X-ray pulsar, recently 
discovered at the center of the LMC supernova remnant 
N157B, close to the 30 Doradus star forming region. 
N157B belongs to the class of the so called Crab-like 
supernova remnants, or plerions, characterized by non- 
thermal spectra and a centrally-filled radio/X-ray mor- 
phology probably due to the presence of a synchrotron 
nebula powered by the relativistic wind from a young, en- 
ergetic, pulsar. Apart from the Crab Nebula, so far only 
three other plerions were known to host pulsars, namely: 
SNR0540-69 (PSRB0540-69), also belonging to 30 Do- 
radus complex and located 15' from N157B, MSH— 15- 
52 (PSRB1509-58) and Gil. 2-0. 3 (PSR J1811 — 1926). 
Thus, with the detection of PSR J0537-6910, N157B rep- 
resents the fifth case of a plerion/pulsar association. 
Pulsed X-ray emission at 16 ms was serendipitously dis- 
covered during a RXTE/PCA observation towards 30 Do- 
Send offprint requests to: R.P. Mignani (rmignani@eso.org) 

* Based on observations collected at the European Southern 
Observatory, La Silla, Chile 



radus (Marshall et al. 1998). Soon after, the pulsation 
was detected in archived 1993 ASCA/GIS data and ad- 
ditional confirmations came from BeppoSAX (Cusumano 
et al. 1998). PSR J0537-6910 takes over the Crab (33 ms) 
as the fastest "classical" (i.e. not spun up by matter ac- 
cretion from a companion star) pulsar. 
The pulsar was identified in the ASCA/GIS with a X-ray 
source detected at the center of N157B and resolved in a 
point-like component (the pulsar and, probably, its asso- 
ciated synchrotron nebula) plus an elongated feature, the 
origin of which is still uncertain. In both RXTE and ASCA 
data, the pulse profile appears characterized by a sharp 
(~ 1.7 ms FWHM) symmetric peak, which shows no ob- 
vious evolution during the time interval between ASCA 
and RXTE observations (3.5 yrs). The period derivative 
of the pulsar (P ~ 5 x 10~ 14 s s _1 ), obtained from the 
comparison of multi-epoch timing (Marshall et al. 1998; 
Cusumano et al. 1998), gives a spindown age of « 5 000 
yrs), similar to the age of the remnant estimated by Wang 
& Gotthelf (1998a), a magnetic field of w 10 12 G, typi- 
cal for a pulsar this young, and a rotational energy loss 
E ~ 4.8 x 10 38 ergs s^ 1 . Substantially the same results 
were obtained from the timing analysis of ROSAT/HRI 
data (Wang & Gotthelf 1998b). 

In radio, PSR J0537— 6910 has been observed between 
June and August 1998 using the 64m radio telescope in 
Parkes but it has not been detected down to an upper 
limit of F14.GHZ ~ 0.04 mJy (Crawford et al. 1998). Al- 
though not really compelling, the present upper limit sug- 
gests that PSR J0537— 6910 is weaker in radio than both 
the Crab pulsar and PSRB0540-69. 



2. Optical observations 

While in radio PSR J0537— 6910 is an elusive target, in 
the optical domain the situation appears more promising. 
Up to now, three of the five pulsars younger than 10 000 
yrs have been certainly identified in the optical (Mignani 
1998), where they channel through magnetospheric emis- 
sion ~ 10~ 5 — 10~ 6 of their rotational energy output. Since 
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PSR J0537— 6910 is very young (~ 5 000 yrs) and, with the 
Crab, it has the highest E, it is natural to assume that 
also in this case a significant amount of the rotating power 
be radiated in the optical. However, given the uncertain 
dependancc of the optical luminosity vs. the pulsar pa- 
rameters, it is difficult to make a prediction on the actual 
magnitude of PSR J0537— 6910. A possible estimate can 
be obtained by a straight scaling of the Pacini's relation 
(see e.g. Pacini & Salvati 1987) i.e. neglecting the depen- 
dance of the pulsar luminosity on its unknown optical duty 
cycle. This would yield V ~ 24.6, after correcting for the 
interstellar absorption Ay ~ 1.3, estimated applying the 
relation of Fitzpatrick (1986) with an Nr <~ 10 22 cm~ 2 , 
measured by the X-ray spectral fittings (Wang & Gotthclf 
1998a). However, we note that the other young (~ 2 000 
yrs) LMC pulsar, PSRB0540-69, with a factor 3 smaller 
E, has a magnitude V = 22.4 with an Ay ~ 0.6 (Caraveo 
et al. 1992). 

Although PSRJ0537-6910 is still undetectable in ra- 
dio, its detection in ROSAT/HRI data (Wang & Got- 
thclf 1998b) reduces significantly its position uncertainty 
down to ± 3"and prompts the search for its optical coun- 
terpart. The scientific case appears similar to the one 
of PSRB0540— 69, also discovered as a pulsating X-ray 
source (Seward et al. 1984), also embedded in a supernova 
remnant (SNR0540— 69) and tentatively identified in the 
optical without the aid of a reference radio position (Car- 
aveo et al. 1992). 

In the following, we describe the results of the first deep 
imaging of the field of PSR J0537— 6910, performed with 
the ESO/NTT. 

2.1. The data set 

The field of PSR J0537-6910 has been observed in three 
different runs between September and November 1998 
from the European Southern Observatory (La Silla). The 
observations have been performed in visitor mode with 
the NTT, equipped with the second generation of the SU- 
perb Seeing Imager camera (SUSI2). The camera is a CCD 
with a field of view of 5.'5 x 5.'5, split in two chips, and a 
projected 2x2 binned pixel size of 0716. The two CCD 
chips are physically separated by a gap ~ 100 pixels in 
size, corresponding to an effective sky masking of ~ 8". 
Images were obtained in different wide-band filters 
(B,V,I) and in the narrow-band H a , with the available 
data set summarized in Table 1. 

After the basic reduction steps (bias subtraction, flat- 
field correction, etc.), single exposures have been com- 
bined and cleaned from cosmic ray hits by frame-to-frame 
comparison. The frames taken through the same filter have 
been registered with respect to each other and combined 
through a median filter. The conversion from instrumental 
magnitudes to the Johnson standard system was obtained 
using a set of primary calibrators from Landolt fields ob- 
served at different airmasses during each night. The formal 



errors in the zero point of the calibration curves are 0.04 
magnitudes in V and B, and 0.03 in I. 
Astrometry on the field has been computed using as a ref- 
erence the coordinates and positions of a set of stars ex- 
tracted from the USNO catalogue. Then the sky-to-pixel 
coordinate transformation has been computed using the 
ASTROM software (Wallace 1990), yielding a final accu- 
racy of 074 on the astrometric fit. 

Fig. 1 shows a 1200s exposure H a image of the 30 Doradus 
region, obtained through the SUSI2 camera. The solid 
square (20" x 20"), located close to the maximum of the 
emission in the H a band just at the center of the star form- 
ing region, includes the X-ray position of PSR J0537— 6910 
(Wang & Gotthelf 1998b). A zoomed I-band image of this 
area is shown in Fig. 2 in negative greyscale. 

2.2. Results 

Few objects are seen close or within the X-ray error cir- 
cle of the pulsar (Fig. 2), including the moderately bright 
(V ~ 19) star #1. However, the crowding of the region, 
together with the irregular sky background conditions, 
prompted us to apply automatic object detection routines 
to search for additional, barely detectable, candidates. 

The object search in the X-ray error circle was thus 
performed using the ROMAFOT package for photome- 
try in crowded fields (Buonanno & Iannicola 1989). The 
ROMAFOT parameters were tuned to achieve in each fil- 
ter a conservative > 5a object detection above the local 
background level. A template PSF was obtained by fit- 
ting the intensity profiles of some of the brightest, unsat- 
urated, isolated stars in the field with a Moffat function, 
plus a numerical map of the residual to better take into 
account the contribution of the stellar wings. To allow for 
an automatic object matching and make the color-color 
analysis faster, all the images have been aligned to a com- 
mon reference frame. As a reference for object detection 
we have used our /-band image, where the effects of the 
local absorption are reduced. The master list of objects 
thus created was then registered on the images taken in 
B and V filters and used as an input for the fitting pro- 
cedure. A carefully check by eye has been performed in 

Table 1. Summary of the available NTT observations 
of the PSR J0537— 6910 field. Columns list the observing 
epochs, the imaging device, the filters, the total exposure 
times in seconds and the average seeing conditions during 
each observation. 



Date 


Detector 


Filter 


T(s) 


Seeing 


Sept 98 


SUSI2 


V 


1800 s 


172 




SUSI2 


I 


2400 s 


170 




SUSI2 


Ha 


1200 s 


1".0 


Oct 98 


SUSI2 


B 


3 x 500 s 


172 


Nov 98 


SUSI2 


V 


2 x 900 s 


078 
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Table 2. Observed, not dereddencd, magnitudes of the 
objects observed inside/close to the X-ray error circle. Ob- 
jects are identified by their labels in Fig. 2. Columns list 
the U-band magnitudes, and colors. Attached photomet- 
ric errors are ~ 0.02 mag and ~ 0.1 mag at V = 18 and 
V = 21, respectively. 



Id 


V 


B-V 


V-I 


1 


18.8 


0.38 


0.45 


2 


21.2 


>1.9 


0.85 


3 


20.4 


0.40 


0.02 


4 


20.7 


>2.4 


0.72 


5 


19.9 


0.32 


0.11 


6 


20.5 


0.58 


-0.70 


7 


20.0 


0.54 


1.58 


8 


19.8 


0.72 


0.62 


9 


19.8 


0.95 


1.11 


10 


20.8 


0.21 


0.74 


11 


20.0 


0.35 


0.36 


12 


17.4 


0.13 


0.27 



order to ensure that all the stellar objects found in the / 
band were successfully fitted in the other images. Apart 
from the ones labelled in Fig. 2, no other candidate opti- 
cal counterpart to the pulsar has been clearly detected by 
our procedure. We just report the possible presence in the 
I-band image of an ~ 22.3 magnitude object (not recog- 
nizable in Fig. 2), right below the detection threshold and 
located nearly at the center of the error circle. However, 
the very low significance of this detection as well as the 
lack of color information prevent us to assess the nature of 
this object and to speculate about a possible association 
with the pulsar. 

The properties of objects #1-#12, i.e. their magnitudes 
and colors (B — V and V — I) are summarized in Table 2. 
According to their colors and brightness, all these objects 
are likely identified as young massive stars. Fig. 3 shows 
the color-magnitude (/ vsV — I) diagram computed for a 
sample of objects selected in a ~ 0'.5x0'.5 surrounding area 
(the region marked by the dashed box in Fig. 1), together 
with the Zero Age Main Sequence track estimated from 
a suitable chemical composition [Z = 0.008, Y = 0.23) 
for the LMC stellar population (Cassisi, private comm.). 
The objects labelled in Fig. 2 and listed in Table 2 have 
been marked by open diamonds. Although broadened by 
the interstellar absorption and shifted redward by the dif- 
ferential reddening, the color-magnitude diagram (CMD) 
of all the stars is indeed consistent with a young stellar 
population main sequence. Thus, the optical counterpart 
of PSRJ0537-6910 is probably too faint to be detected 
against the high background induced by the supernova 
remnant and the absorption of the embedding HII region. 
Using the template PSFs computed from each image, ar- 
tificial stars tests have been run to estimate the V, B, and 
/ magnitude limits of our images. With the flux normaliza- 



tion left as a free parameter, artificial stars have simulated 
and added to the corresponding images at 100 different po- 
sitions randomly selected inside the error circle. Thus, the 
detection algorithm has been run in a loop for the above 
number of trials, with the flux normalization adjusted to 
allow for a 3cr detection in each filter. Averaging over the 
number of trials, we have found 3ct detection limits cor- 
responding (within 0.2 mag) to B ~ 23.2, V ~ 23.4 and 
I ~ 22.4, which we have taken as an indication of the 
limiting magnitudes achievable in each band. 

3. Conclusions 

We have performed deep optical observations to search 
for the optical counterpart of the isolated X-ray pulsar 
PSR J0537— 6910. However, none of the objects detected 
close to/inside the X-ray error circle stands out as a con- 
vincing candidate. The marginal detection of a / ~ 22.3 
object at the center of the error circle must be regarded 
as tentative and is in need of future confirmation. The op- 
tical counterpart of PSR J0537— 6910 is thus unidentified 
down to a 3a limiting magnitude of ~ 23.4 in V. Our re- 
sult is in agreement with the upper limits recently derived 
by Gouiffcs & Ogelman (1999) on the pulsed optical flux. 
At the distance of 47 kpc estimated for the host remnant 
N157B (Gould 1995) and for the assumed interstellar ab- 
sorption (Ay ~ 1.3), our upper limit corresponds to an 
optical luminosity L opt < 1.3 x 10 33 erg s _1 . This implies 
that PSR J0537— 6910 is, at best, of luminosity compara- 
ble to the ones of the Crab and PSRB0540-69, in line 
with the predictions of Pacini's law. Although interesting, 
this upper limit is not stringent enough to put strong con- 
straints on the evolution of non-thermal optical emission 
of young pulsars. Together with the recent upper limit ob- 
tained for PSRB1706-44 (Mignani et al. 1999), the mea- 
surement of the optical luminosity of PSR J0537— 6910 
would be crucial to smoothly join the class of the very 
young (~ 1 000 years) and bright objects with the class of 
older, Vela-like (~ 10 000 years) ones, for which the op- 
tical output is ~ 4 orders of magnitude lower (Mignani 
1998). 

As in the case of PSRB0540-69 (Shearer et al. 1994), 
time-resolved high resolution, imaging, possibly exploit- 
ing the more accurate X-ray position available from future 
Chandra observations, would be the best way to pinpoint 
and identify the optical counterpart of PSR J0537— 6910. 
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Fig. 1. SUSI2 H a image of the 30 Doradus star forming 
region. The location of PSR J0537— 6910 is near the center 
of the small ~ 20" x 20" solid square. The larger (0.'5 x 0.'5) 
dashed box marks the area in which the search of stellar 
objects has been performed (see Sect. 2.2). 
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Fig. 2. ~ 20" x 20" I-band image of the PSR J0537-6910 
surroundigs (North to the top, East to the left). The circle 
(~ 4") marks the uncertainty region associated with the 
pulsar position as resulting from the combination of the 
error (~ 3") on the boresight-corrected ROSAT/HRI co- 
ordinates (Wang & Gotthelf 1998b) and the global error 
budget of our astrometry. Objects detected close/inside 
the error circle are labelled 1-12. 
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Fig. 3. Color-magnitude diagram of the objects found in- 
side the dashed box in Fig. 1 (crosses) and the objects la- 
belled in Fig. 2 (open diamonds). The solid line in the plot 
corresponds to the Zero Age Main Sequence computed for 
stars with a chemical composition Z — 0.008, Y = 0.23 
and masses in the range from m = 1AM® to m = 18M 
(Cassisi, private comm.). 



